function [ vh ] = fmg( fh, current)
%FMG Summary of this function goes here
%   Detailed explanation goes here

global P;
global k;
global v0;

% Step 1: If current iteration is not the last one:
if current < k
    f2h = P{current}'*fh;
    v2h = fmg(f2h, current + 1);
    % Step 2: correct
    vh = P{current}*v2h;
else
    vh = zeros(length(fh), 1);
end

% Step 3:
for i = 1 : v0
    vh = vcycle(vh, fh, current);
end

end

